Decay of Z-string due to the fermion emission. 



A. A. Kozhevnikov 

Laboratory of Theoretical Physics, 
S.L. Sobolev Institute for Mathematics, 
630090, Novosibirsk-90, Russian Federation * 
(February 1, 2008) 



Abstract 



The question of classical topological stability of the gauge vortex defects 
is reanalyzed upon taking into account the quantum perturbative effects. The 
purely Abelian Higgs string remains stable while the coupling of the effectively 
Abelian Z-string (with the fixed zero upper component of the scalar doublet) 
to the fermions of the minimal standard model is shown to result in its decay, 
via the fermionic pair emission. The decay rate is scaled as 10~^ of the 
decay rate of boson. The influence of the surrounding charge-asymmetric 
fermionic matter is considered, demonstrating the dependence of the decay 
rate on the contour shape and its suppression at sufficiently large charge 
asymmetry. 
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I. INTRODUCTION 



The interest to the electroweak (EW) defects in the form of Z-strings [||] that could 
be produced at the phase transition in early universe is attributed, in particular, to their 
possible role in the processes of baryogenesis at the EW scale p[. First, they are treated 
as possible carriers of nonzero baryon number so that the change of this number 

AB is related to the change of the Chern-Simons (helicity) number Ancs of the Z-string 
configuration via the integrated anomaly equation for the baryon number current ||^: 

—2 

= Ancs = ^fj^ cos(2^iy)A j d'xZ • (V x Z), (1.1) 

provided the background electromagnetic field is absent. Here g and g' are, respectively, 
the SU(2) and U(l) coupling constants, g = \J g^ + g''^, 9w is the Weinberg angle, and Nf is 
the number of the fermionic families. The last integral in Eq. ( |1.1| ) is recognized to be the 
helicity number hz of the EW Z-string configuration. Second, EW strings were suggested 
to be the source of nonequilibrium ^ required in any model of baryogenesis, to produce 
the baryon density via the evolution of the EW string network. 

Although the simplest variant of Z-string solution is known to be unstable [jl| towards 
to the development of either the upper component of the Higgs field doublet @,|], or to W 
condensation |T^, these instabilities can be possibly neutralized by particles bound to Z- 



string [mTT| or by imposing external magnetic field [|T2| , respectively. The two Higgs doublet 



generalization of the Minimal Standard Model (MSM) admits classically stable Z-string |]13 

In the papers cited above the possible quantum properties of Z-strings were neglected 
despite the fact that these defects possess the microscopic transverse size of the order of 
~ rriz^ and could behave similar to bosons and, in particular, become unstable via the 
emission of fermions. The latter point should be clarified from the very start. Indeed, the 
Z-string with the fixed zero upper component of the Higgs field doublet can be treated as 



an effectively Abelian Abrikosov-Nielsen-Olesen (ANO) string [14| in the Higgs model with 
the action 



S= (Tx 



Fl + \{d, + igA,) 




(1.2) 



This string is known to be topologically stable. The formal proof is reduced to the statement 
that the homotopy /(p)[(l — M)exp(i0) + u] interpolating between the two static configu- 
rations of the scalar field (f) with different winding numbers n = 1 at m = and n = at 
M = 1, is not admissible since the energy barrier A_E oc v?{l — uY dppf^{p), separating 
these two configurations with the finite energy per unit length is infinite. Here g is the 
U(l) coupling constant, the profile of the scalar field /(p) has the asymptotic / 
p — s> oo, 6* is the polar angle in the xy plane, p = ^/x^ + y^. In quantum case one examines 
the path integral 

r r<f>(tf)=4>f 
{^f\^,) = / d[ct>f]mct>f] / d[ct>\eM^S[m^[Hdm. (1-3) 

{(f) is the shorthand notation for all relevant fields) which determines the overlap of two wave 
functionals. The latter is claimed to be vanishing, if the fields (pi and belong to distinct 
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homotopical classes |]T5[ characterized in the present case by different winding numbers of 
the scalar field, since the infinite energy barrier makes the exponent to be rapidly oscillating. 
How is evaded this topological veto in the course of particle emission? 

To this end one should emphasize the following. First, the above homotopy determines 
the field configuration that does not obey the equations of motion, except for the border 
values u = and u = 1, provided u is taken as dynamical variable. Admitting it as 
the quantum trajectory which does not oblige to obey the classical equations of motion, 
one encounters the infinite energy barrier making the amplitude vanishing. So, strictly 
speaking, only this specific type of the field deformation is not allowed quantum mechanically. 
On the other hand, the time dependent winding number n{t) could interpolate between 
different topologies without going beyond the class of fields with finite energy per unit 
length, = T] exTp[in(t)9]/ ^/2. This is still the solution of the equations of motion, provided 
the time component of the vector potential is induced via the condition of the vanishing 
covariant time derivative of the scalar field. (See Sec. ^ below). Classically, however, the 
winding number being the degree of mapping should be an integer which would seem to leave 
only the known possibility of the conserved n. One of the purposes of the present paper is to 
point out that this conclusion can be evaded by treating this number as quantum dynamical 
variable n = n{t). To this end we go beyond the approximation of the classical background 
field and quantize the vortex field configuration. The latter, in the approximation of the 
Higgs boson mass being much greater than the gauge boson mass called sometimes as the 
London limit, is known to be characterized solely by the spacetime dependent phase of the 
scalar field. The demand of the scalar field to be single valued should be imposed then only 
on the classical field configuration. Precisely, the expectation value of n in the quantum 
state \1/ should be an integer. The wave function of this state \1/ will be found explicitly. 
The next step is to study the evolution of this wave function under the influence of the 
interaction of the vortex background with the propagating excitations. It will be shown 
that the wave function in its dependence on time and the winding number acquires the 
components peaked at successively diminishing winding numbers which is interpreted as 
violation of the topological stability of the classical string configuration. The lifetime of the 
string with the unit winding number turns out to be proportional to the lifetime of the gauge 
boson of the model, with the factor depending on the string contour. The purely Abelian 
string in the model ( |1.2| ) survives this mechanism due to kinematical reason, since the net 
mass of the final particles entering the vertex of their interaction with the string background 
is greater than the energy splitting of the quantum levels between which the transition 
occurs followed by the particle emission. The emergence of the levels results from treating 
the winding number as quantum variable [|T^]. It should be emphasized that this type of 
instability is by no means include the transition via the infinite energy barrier mentioned 
above. 

In order to fix the notations, let us consider the neutral current piece of the lagrangian 
density of MSM assuming for a while a single fermionic family: 

•^NC -^boson ~l~ -^fermion ~^ -^Yukawa; 

where the bosonic part is 

^boson = -\zl + \{d, + '-gZ,)<!>\' - ^-X'm' - \r,')\ (1.4) 
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and the upper component of the Higgs doublet is taken to be zero; the fermionic part is 



+ 



iil)+R d -ig{-Q+R^)Z tp+R + iip_R d - ig{-Q-R^)Z tp.R 



while 



:i.5) 



Yukawa 



-h^{ip_Lip-R^ + c.c.) - h+{i)+L'ip+R^* + c.c.) 



;i.6) 



gives the mass to the fermions and describes their interaction with the Higgs boson. Note 
that /i+ = in the case of leptons. Hereafter d = 7^9^ etc., Z^^^, = d^^Z^ — d^Z^, ^ = sin^ 9w, 
T3 is the third component of the weak SU(2) isospin, and Ql = diag^Q^i, Q-l) is the charge 
matrix of the left fermions, Q±l(r) being corresponding electric charge. 

The subsequent material is organized as follows. Section y contains the derivation of 
the effective reduced action for the gauge vortex state in terms of the winding number, 
including the discussion of the contour motion omitted in Ref. []TB[. It is argued there why 



the dynamics of the radial part (modulus) of the scalar field is inessential for evaluation of 
the lifetime of the string in the London limit. The quantum state of the string configuration 
and the criterion of its stability against the particle emission is discussed in Sec. |T| . The 
purely Abelian Higgs model admits the string stable against this mechanism. This is not 
the case for effectively Abelian Z-string which possesses the coupling with almost massless 
fermions. The calculation of the corresponding decay rate of the state of Z-string having 
arbitrary geometric shape, including the parity-odd contours, is presented in Sec. with 
the taking into account of the effects of the nonzero fermionic density. The damping of the 
energy, Z-flux and helicity (Chern- Simons) number of Z-string is discussed in Sec. |V[ Sec. 
V| is devoted to the discussion of the vahdity of assumptions adopted in the paper, and to 
the conclusion drawn from the present study. 



II. THE ACTION OF Z-STRING IN TERMS OF THE WINDING NUMBER 

The fact that the Z-string solution is an embedding into SU(2)xU(l) gauge group 
of the known U(l) ANO vortex solution [compare Eqs. ( |1.2D and ( |1.4D ] permits one to 



refer to our earlier derivation |jT6| of the reduced effective action. When so doing, it is very 
comfortable to use the space Fourier transforms of the fields, which are elementary functions 
rather than the special ones used in writing down the original solution 

The nonstationary field configuration of Z-string is expressed through the spacetime 
dependent phase x = x(^)^) oi the Higgs field $(x, t) = r]exp{ix)/ in the London limit. 
This is the limit of rriH ^ rriz, InmH/mz is also large, where mz = gvf^ ^ind mn = A?7 are 
the masses of the Z° and Higgs bosons; A and ///a/2 are the Higgs field self-coupling and 
magnitude. Let us remind briefly why the dynamics of the phase x of scalar field, not 
the dynamics of its radial part (modulus) /, is the only one essential in this limit. To this 



end one should rewrite the action Eq. ( |L2|) in terms of these variables, 

S = ld''x [^i-dtA - VA,y - i(V X A)2 + f[{x + gAtf - (Vx - gAf 
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+/^-(V/r-y (101^-^1 }, (2.1) 

and to convince that the contribution of the radial part to the path integral Eq. ( p. .31 ) 
is factored out in the London limit. Indeed, an exact factorization could be broken, in 
principle, in the following two cases. First, there is the mixed third term in the first line of 
Eq. But for large distances p ^ mjj^ one may set / ~ i]/V^, which leaves only the 

phase, while at the short distances p < mjl the amplitude of the Higgs condensate behaves 



as f ^ p^ IjTjl , so this term gives negligible contribution. The conclusion remain valid in the 
case of oscillating n (see below), since one should average the contribution over the period 
of oscillations after which it is not the winding number itself but its amplitude that enters 
into expression leaving the suppression at p ^ intact. Second, the term is proportional 
to at the short distances, however, its contribution to the action is suppressed as 
as compared to the logarithmically enhanced (In ^ ^ 1) contribution of large distances 
p ^ rn'^ coming from the terms containing the phase. So, one can ignore the details of the 
Higgs field profile /, taking it to be uniform r]/\/2 in all coordinate space except the vortex 
line where it approaches zero at characteristic distances ~ ■ The contribution of the 
radial part is then factored out as redefinition of the measure in the path integration over 
nit). 

The equation for the Z-magnetic field, 

2rn? 

V X Hz = — ^ Vx - m|Z, (2.2) 
is solved to give the Z-magnetic field strength, 

Hz(k,t) = ^ ■ J^} , /daXlexp(-zk ■ X,), (2.3) 

and the vector potential: 

Ann / 1 1 \ /■ 

Z(k, t) = —^-[^,- ) f ^^^[k X XI] exp(-zk ■ X.). (2.4) 



Here the integral over a comes from the equation for the phase x read off from Ref. |]T7 
with the proper continuation to the Minkowski spacetime: 

V X Vx(x, t) = 2nna j> daX^^^^) _ x,(a, t)], (2.5) 

where Xq = Xa(cr, t) is the evolving closed string contour a parametrized by the arclength 
a. Hereafter the prime over X will denote the derivative with respect to corresponding 
parameter along the contour, while the overdot will do the time derivative. The case of 
many contours is embraced by taking the sum over individual contributions in the right 



hand side of Eq. (|2.5| ). It is argued in Ref. |]T8l that possible fermionic zero modes [|l^] do 



not perturb the string profile. Recall that the winding number = na{t) of the scalar field 
is directly related with the number of quanta of the magnetic-like flux, in the present case 
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the Z-flux, via the condition of the vanishing covariant derivative of the Higgs field deep 
inside in the Higgs condensate, 

Z-dl=-fvx-dl = (^ona, (2.6) 

where 0o = y is the quantum of Z-fiux. It is essential that the first term in the parentheses 
of Eq. ( |2.4| ) is in fact a pure gauge one. Indeed, the Fourier component of 

2 

v(x,t) = -Vx 
9 

found from Eq. (p.5|), is 



v(k, t) = -^j dai[\<. X XI] exp(-zk ■ XJ. (2.7) 

Locally, in a plane transverse to the tangent vector X^, Eq. ([^.T]), after going back to the 
coordinate space, leads to the transverse components of v to be 

fi(x, t) oc Sij^dj In IpI = c^iatanpa/pi; X^ ■ p = 0. 

To specify the dynamical part of the problem, one should write down the Z-electric field 
Ez = -VZt - 9iZ, where 

Zt = --dtx (2.8) 
9 

replaces the condition Zt = appropriate in the static case. One has 

Ez = -Vdtx - dtZ = -{Vdt - dtV)x + dt i-Vx - Z 1 . (2.9) 
9 9 \9 J 

The commutator of the derivatives is nonzero in view of the singular character of the phase 
X [17], so the Fourier component of E^ becomes 

An 
9 



An r 

Ez(k, t) = —Ua f da{Xa x X^) exp[-zk ■ X,(a, t)] 

U J 



Note that the vector potential Z and the magnetic type field strength can also be 
expressed through the gradient of the singular phase x 

Z(k,t)= fl- j^^^l v(k,t). 



m 



i^z{k,t) = -^^[kx^r{k,t)]. (2.11) 

In what follows we will omit the encounters of the nearby string segments. Important are 
in the processes of the string rearrangements, they cannot be described in the framework of 
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the London limit and demands the numerical integrations of the full set of the equations of 
motion. Substituting Eqs. ( p.3|) , (|2.4| ) and (|2.10|) into the lagrangian Lboson = / '^^a^'^boson 
one obtains, with the help of the relation 



J rf^xHKx) = J rf3A;|H2(k)|V(27r) 



etc, the expression for the action of the single gauge vortex: 

'S'vortex = 2(^2tt)^ J (k^ + TTi^ ^ J J ^^i'^^2 exp{ik ■ [K{ai) - X(o-2)]} 
x{[nie - m|(k2 + m|)n^][X'(ai) ■ ^'(a^)] 
+nl{k' + 2m|)(k ■ [X(ai) x X'(ai)]) ■ (k ■ [X{a2) x X'ia^)]) 
+n>|[X(ai) X X'(ai)][X(a2) x X'ia,)]}. (2.12) 



Let us make a step apart and show with the method similar to those of P. Orland, Ref. |17 



and pOf how the known Nambu-Goto (NG) action results from Eq. (|2.12| ). To this end one 
should take the limit of the fixed winding number and set the Z-boson mass mz — > oo 
before the momentum integration. Then the term oc drops out, and the action becomes, 
in the gauge X° = t = t, 



k2 



^NG = f / dh^dh^S^'^Xis,) - X{S2)]{-X'{S,) ■ X'{S2) 

+ [X(si) X X'(si)] ■ [X{s2) X X'(s2)]}, (2.13) 

where Si^2 = sf^2 — (''"1,2! 0^1,2) is the two- dimensional vector. Using the Gaussian regular- 
ization of the 6 function and the expansion 

X(s2) - X(si) + {S2 - Si)^9^X (2.14) 

20] valid under the condition |X"(o")| ^ m^, one obtains 



^NG = \ / dh,d^zeM-i^2''^z''dAX^dBX,){-X!' + [X X Xf ) 

j dh^detdAXi'dBX^, (2.15) 



47rA2 

where ^ 00 is an ultraviolet cutoff, Oa = d/dz^ and dei^AX^'^BX^ = -X'^+[X x Xf 
in the chosen gauge. Up to an overall factor, the last equality in Eq. ( |2.15| ) is recognized to 
be the NG-action. 

Coming back to the case of large but finite mz, one obtains the action as the sum 
Ea Cortex, where 



42tex = ^\n^Jdtfda{inl- m|n^)Xl^ + m|n^[X. x X'f/cl} . (2.16) 



Here the first and the third terms come from the electric-type field, while the second term 
comes from the kinetic energy of the scalar field, thus demonstrating an intimate interplay 
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of the vector and scalar fields in the string solution. The energy of Z-magnetic field is not 
enhanced logarithmically in the London limit and by this reason is dropped. Further, 

2 Ami- , tuh 
cl= \ ln^<l 

is the velocity squared which characterizes the classical string motion in the case of finite 
masses, the factor 0(1) reflects the ignorance of the true Higgs field profile, and rriH appears 
as the natural upper limit of the integration over momentum. We retain only the terms in 
the action that refer to the string background. Omitted are the terms which correspond to 
the propagating Z and Higgs bosons. They result in the renormalization of the parameters of 
the action Eq. ( p.l6|) and do not contribute, due to the energy conservation, to its imaginary 



part (see below), which will be further of our main concern. Eq. ( p.l6| ) is valid in the case of 
large but finite m^, so the terms corresponding to the interaction are exponentially small, 
and only nearby segments of the string contour give an appreciable contribution to the 
integral over arclength. 



III. QUANTUM STATE OF THE STRING CONFIGURATION AND THE 
CRITERION OF ITS STABILITY AGAINST THE PARTICLE EMISSION 



It is seen that under the condition |Xa| <C Cq assumed hereafter, the dynamics of the 
winding number becomes decoupled from the contour dynamics and is governed by the 
oscillator-like action, 



(a) 



vortex 



2 2\ 



where the frequency and effective mass, in the gauge X'^ 



SnL iriH 
m . 



mz 



(3.1) 

1, are, respectively, oj = mz and 

(3.2) 



Hereafter L = § da is the length of Z-string in the chosen gauge. The quantization subjected 
to the constraint {il)\n\ip) = tiq, where no is an integer, is performed with the help of the 
indefinite Lagrange multipliers and gives the energy levels of the system 



EN = e,Lnl + mziN+-) 



(3.3) 



where = Airm'^g ^ In mn / mz is the energy per unit length of the vortex with the unit 
winding number, = 0, 1 ■ ■ ■ (do not confuse with the winding number n) pT|. The wave 
functions are the oscillatory ones, tpNi^) = ip^^'^^^' 



pgyLy/^ 1 



exp 



[n - no) 



mz 



■{n - rioY 



H 



N 



mz 



(3.4) 



(i^AT is the Hermit polynomial) displaced to no. They are sharply peaked at the integer 
numbers. Specifically, the wave function of the vacuum state without string is ip, 
irrespective of the string contour X and the radial part of the scalar field. 



(osc) 
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Now one can evaluate the overlap integral (|1.3| ) for the string with different winding 
numbers, rii at the moment t = and Uf at some later moment t > 0. The string will be 
assumed to be in the ground states, = 0, so that its wave functions at those moments 
are ipi = iljQ'^'^\n — Ui) and '?/'/ = ipQ'^^\n — Uf), respectively. The result depends crucially 
on the spectrum of the propagating excitations, since the string interaction with the latter 
may result in an imaginary correction to the energy levels Eq. ( p.3|) . As it will be shown 
below, such a correction emerges as the imaginary part of the frequency of the ra-oscillator, 
rrizN —>■ [mz — iV /2)N . Using the known expression for the propagator of the damped 
oscillator, one obtains the probability amplitude of the transition with the change of the 
winding number: 



(f\i) = exp I — nf + nl — 2ninf exp(— rt/2 — imzt) \ . (3.5) 

I 2mz L J J 



The expression for the probability averaged over the period T = 27rm^^ looks as 

{wfi)=Afexp —^^(n^ + nl) Iq ^^-^nj^f exp (—rt/2) 
L niz J L mz 

~ AT exp [nf - niexp(-n/2)]'| (3.6) 

[lo{z) is the modified Bessel function of order zero], where the contribution of the path 
integral over the radial part of the scalar field, as is argued in Sec. |T|, is taken into account 
in the form of the constant normalization factor A/". 

If the damping were absent, F = 0, the expression for the relative probability would 
have a sharp maximum at rij = nf, thus reproducing at the quantum level the topological 
conservation of the string winding number. On the contrary, at T 7^ the maximum 
is achieved at the configurations with the decaying winding number. The character of 
the decay is established by reconciling the integers Ui and nf with the continuous time 
dependence in Eq. ( |3.6| ). Indeed, somewhere within the time interval At from the start, 
such that -2T^^ ln(l - 1/2^^) < At < ~2T-^ ln(l - 3/2^^), the nearest inte ger nf for which 
the probability is nonzero, is rij = rij — 1, and so forth down to ni = 1. In particular, the 
lifetime of the U(l) gauge string with n = 1 is 2r~^. The classical topological veto is evaded, 
because due to the damping the wave functional of the initial string configuration acquires 
the components of the vacuum state. Note that the topology change is faster for initially 
multiple winding numbers, since the time of awaiting of = rij — 1 is At ~ l/Erij ^ 
p^ , so that the total decay time is T ~ 2r"^ SfcLi ~ 2r~^ Inrij. The last approximate 
expression has the same form as if the winding number would change continuously in the 
course of the decay. 

Let us confront the result of Eqs. (|3.5|) and (|3l6|) with the usual approach to the string 
stability. What is implied in the present approach is the preparation of the string in the 
quantum state in which the expectation value of the winding number is kept fixed. It 
is perfectly the case that takes place for the Abrikosov string in type II superconductors 
immersed into external magnetic field, while no specific mechanism has not yet been devised 
for the Z-string. Now remove the constrain sustaining the string. Then the state with the 
fixed winding number becomes the superposition of the stationary states of an unconstrained 
oscillator, so that the amplitude of finding it in the state with different expectation value nf 
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of the winding number at a later moment is given by Eq. ( p.5|) . The feature of the present 
approach is thus the quantum mechanical treatment of the background field configuration. 
It should be recalled in this respect that the notion of the classical background field is 
itself an approximation aimed to represent the condensate of indefinite number of quanta 
of corresponding quantum field. One of the possible representatives of such a condensate 
is the coherent state, in which the expectation value of the field operator, in the present 
case the phase of the Higgs field, is nonzero. The wave function ipl)'^'^\n — no) used in 
the above calculation is just the wave function of the coherent state in which the winding 
number of the scalar field is the integer uq. In this sense the usual topological stability of 
the classical configuration would correspond to the permanently constrained string resulting 
in the amplitude {f\i) oc 5„.ri^ irrespective of the damping rate T. The role of the damping 
in this case is reduced to relaxing the field configuration to that determined by external 
conditions. 

To establish if F is zero or not and hence the condition of the string stability, one 
should examine the spectrum of the propagating excitations and their couplings with the 
string background. Expanding the action of purely Abelian Higgs model Eq. ( |1.2|) into the 
contributions of the string background and the excitations, which are the massive neutral 
Higgs (pn and vector A bosons, one can show that the only relevant coupling is the string—*- 
ifH + A one. Since the emission of the excitations occurs in the course of the transitions 
between the levels spaced, with the evident replacement, by ~ m^, its contribution to F 
is forbidden by the energy conservation. Thus, the stability of the Abelian Higgs string in 
the model Eq. (|1.2|) with the neutral scalar field is merely the kinematical consequence of 



rather limited spectrum of excitations It well may be not the case in other models. 

In particular, the effectively Abelian Z-string with the fixed upper component of the Higgs 
doublet illustrates this since possesses the interaction with almost massless fermions. 



IV. EVALUATION OF THE DECAY WIDTH OF Z-STRING. 

The fermionic loop correction to the action Eq. (|3.1| ) can, in principle, be evaluated by 
the integrating out the propagating fermions from the total action of the standard model. 
However, we need only the imaginary part of the resulting effective action. Since this 
imaginary part arises due to real intermediate states on the mass shell, one can use the 
unitarity relation for its evaluation and, in turn, Imu, allowing for the quantum transitions 
between the energy levels of the oscillator Eq. (|3.1|). 

In order to obtain the lagrangian of interaction of Z-string with the physical fermions, 
one should rotate away the phase x from the mass term of the Lagrangian of the standard 
model Eq. ( |1.6|) , with zero upper component of the Higgs field. It can be accomplished by 
the phase rotation of the chiral fermions, 

^±R ^±R exp (2ix<5±i?0 5 (4-1) 
with the phases proportional to their respective Z charges, 

Qz = T3- Qsin^ 6w] 
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Q and T3 standing, respectively, for the electric charge and the third component of weak 
isospin of the chiral fermion. The pure gauge term of the background field Z presented as 
the first term in the parentheses of the expression for Z in Eq. absorbs the gradient 



of the rotation phase thus leading to no shift of the lagrangian. The term describing the 
interaction of the single chiral fermion, say the left one, with the string becomes 

>Cint = gi'LZ'-^^Qz^L, 

where the Fourier transform of the space components of the short range piece of the back- 
ground Z-field Z^^-* is given by the second term in parentheses of expression for Z in Eq. 
( p. lip , while the time component is zero in the London limit. Note that taking the limit 
rriz —>■ 00 results in the decoupling of Z-flux from external fermions. 

The matrix element of the emission of a pair of chiral fermions by Z-string in the course 
of the quantum transition between the levels and — 1 shifted by mz is 

Mn = -2ngQz5imz - si - £2)ml(pi)7 • ^vl{p2){N - l|n,|iV), (4.2) 

where the above Z*-^-* is factored out as naV(pi + P2), and the recoil is evidently negligible. 
Here 7 is the spatial Dirac matrix, pi = {si, pi) is the four-momentum of the final fermion 
i = 1,2, and 

{N - l\na\N) = {N/2M,mzf'\ (4.3) 



as is evident from the oscillator-like character of the action Eq. (|3.1|) . The following cal- 
culation of the width A^v of the Nth level includes the effects of the particle-antiparticle 
asymmetry in the surrounding fermionic matter by means of the Fermi blocking factors in 
the expression for the density of final states. Inserting the identity 1 = / d^k5(}<i — pi — P2), 
one obtains 



= 7T~^ / ^ ^ ^'^^z ~ ^1 ~ ^2)5(k - Pi - P2) 

8M^mz[27!-)^ J 6162 



X 



E l^dPih ■ ^VLiP2)ni - f{ei - /x)][l - f{e2 + /i)] 



spins 



loirM^mz J (37r)'^|k| J(mz-|k|)/2 
+ (k2 - ml){2e - mz)tk ■ [V x V^J/k^} [1 - /{s - - f{mz - e + ^)]. (4.4) 

Here T, /i are, respectively, the temperature and chemical potential, V = V(k), /(cu) = 
(exp ^-|- 1)~^ is the Fermi distribution, the fermions are taken massless, and the temperature 
dependence of niz is assumed to be included. The linear dependence of A^r = A^F on 
means a nonzero imaginary part of the frequency Xmio = — F/2, hence the damping of the 
quantum mechanical probability, with the rate F. 



Internal integral over e in Eq. ( [4.4| ) can be written as 

»(»Tiz + |k|)/2 
/{mz-|k|)/2 



G(/x,|k|)= / d£---=|Vp(k2 + m|)/i(|k|) 

i{mz-|k|)/2 

+2{2mz\V\^ + (k^ - ml)tk ■ [V x V*]/k^}l2{\k\), (4.5) 
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where 



im)- ^ In ^hf + ch^ 

1 - eM-mz/T) chf + ch^ ' 



/2(|k|) = -sh£exp^ / de-, ^ ^. (4.6) 

211 \) 2 T ^2Ti-|k|)/2 (ch^ + ch^) (ch^ + ch^) ^ ^ 

Because of V is rather comphcated function of the momentum, further evaluation of the 
decay width cannot be performed exphcitly. Instead, we will obtain useful approximate 
results valid in the situations of the physical interest. 

i) Empty space, yU = 0, T = 0. The charge conjugation (C)-even, parity (P)-odd, CP- 
odd structure ik- [V x V*]/k^ drops after the integration over final states as it should, since 
no CP nonconserving effects remain in the charge symmetric situation, in the absence of an 
explicit CP nonconservation in the lagrangian. The necessary expression for |V(k)p, 

\y(k)\'=(—] /rfMa2(x^x-,)expHk-(Xi-X2)], 

where Xi = X((Ti) etc, can be obtained for sufficiently smooth contours whose curvature 
satisfy the condition | X" | / mz <^ 1 • Then one can use the expansion 

X(a2) = X(aO + ^X'(aO + ^X"(ai) + |^X'"(ai) ■ ■ • , (4.7) 

where z = o"2 — ai, to show that upon neglecting the terms with the second and higher 
derivatives of the contour the following approximate expression holds: 




|V(k)r^2vr(^) (k^^^/^^^(k-X'.). (4.8) 



When obtaining Eq. ( [4. 8] ), the integration over z can be extended to ±oo. Choosing the local 
(at given a) coordinate system k = (k_L, k ■ X'), the k integration is easily performed. The 
summation over all fermionic species, with the expression for the total width of Z°-boson 



gives the damping rate 



Hereafter the notation 



rz = frnz E Ql/24vr, 

fermions 



r = fi:(l-ln2)rz. (4.9) 



Alnmn /raz 

is used. The length of the string L drops from the final expression, since the factor L^^ 
coming from the probability of the quantum jump |(A^ — l|na|A^)p [see Eqs. ( |3.2| ) and 



LSI)], is cancelled by another factor L [see Eq. ( |4.8| )] arising due to incoherent emission of 
a fermionic pair by arbitrary point along the string. Since In m^^/m^ is conjectured to be 
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large, say 3, the damping becomes about 10~^ of the decay rate of Z° boson. Note that the 
emission of the Higgs boson and of a pair of Z-bosons is forbidden by the energy conservation, 
since in the lowest order adopted here the energy released in the decay is m^. The photon 
emission occurs only in higher orders, hence corresponding partial rate is suppressed. 

ii) Cold, strongly charge-asymmetric fermionic matter, |yu|/T ^ 1, > uiz- This is the 
case of the strong degeneracy of the fermions, so one can write 

,1 _ - _ e.p,™;!'^':!;,';?: ;r V, < o. ^^-^^^ 

Consider for the definiteness the case /i > 0. Then the expression for G(/i, |k|) reads 



G(^,|k|)^Texp!^M^||Vp 



\2 |k| I ^2 -lkl 



|k| + m^) exp — — (|k| — m^) exp 



+2|k|sh^ ■ (k^ - m|)zk • [V X V*]/k2| . (4.11) 

Since the dominant contribution comes from |k| ~ niz, and ruz ^ T, one may keep only 
the rising exponents in the calculation. Furthermore, all the polynomial |k|-dependent 
expressions in the denominators of the V-dependent expressions can be approximated by 
their values at |k| = rriz- The evaluation of the CP-odd contribution uses the expression 



^k-[vxv*] f^TcY j:^ zk-[x;xx'2] , ^ , 

' ^daif da2 (]^2^^2)2 exp(-«k ■ X12) 




where X12 = X(cri) — X(cr2), together with the dropping of the terms, which are suppressed 
by the powers of m2|Xi2| ^ 1. The integration over k and the inclusion of the case /i < 
results in the damping rate: 

T^J^)\eM^^))U + ^)T,, (4.12) 

KnizJ T \ mz J 



where 

((■■■)) = E(---)Ql/EQ 

/ / 

is the average over Z-charges, and 

X12 ■ [X']^ X X2 



2 

Zi 



Fi = — f daida2 sinm^|Xi2| (4-13) 



ttLJ X 



is the parity-odd factor which depends on the geometry of the string contour. 

The damping is completely prohibited in this case. In fact, the energy of emitted fermion 
spreads from to mz, so at > niz it necessarily hits the filled states inaccessible by 
the Pauli exclusion principle. The conditions of the complete suppression are not likely 
to be achieved in the universe at large scale. To be sure, the chemical potentials of all 
fermions should obey the condition > mz- But only for neutrinos there is a rather high 
cosmological upper bound [I]i/(/^i^/^)^]^''^ < 45 for the net contribution and considerably 
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more stringent bound from nucleosynthesis, \^yJ\/T < 0.2, for the electronic neutrino p5 
The baryonic and leptonic asymmetries are small, ~ 10~* — 10~^°. So the averaged fermionic 
asymmetry of the universe is insufficient to prohibit the decay strongly. Yet the smaller scale 
inhomogeneities in chemical potential of all types of fermions at the EW epoch could not be 
excluded by this argument. 

iii) Cold, weakly charge-asymmetric fermionic matter, \fi\ ^ T ^ mz- One has 

G{fx, |k|) ~ |k| || V|2 (k^ + m| + Sfxmz exp 

+4/xexp ^ ■ ch^ ■ (k^ - ml)^k ■ [V x V^J/k^j , (4.14) 
which results in the damping rate: 

with the same geometric factor Fi as in previous case. 

iv) Hot, weakly charge-asymmetric fermionic matter, n/T -C 1, mz/T <^ 1. Here the 
above G is 



G(^.|k|)=.l|k|{|Vp[k=(l + ?g^^U 



rrir 



12T2 



(k'-m|)ik- [Vx V*]/kn, (4.16) 



The evaluation of the CP-even contribution looks similar to the previous cases, while the 
CP-odd term requires some care. One has 



/g— |k|)^ . ^ (f / ..,(|x; X x; 



(27r)3|k| 12T2 \^ g 



2\ ■ '^12 



x(-i) /^e(m2-|k|)|i-^i^+ ^'"'^ 



X exp(— zk ■ X12) 



k^ -|- m| (k^ -|- m|)2 



When evaluating the integral over m^-dependent terms, one should have in mind that the 
upper integration limit can be set to infinity. Indeed, the oscillating factor sin|k||Xi2| 
resulting from the integration over the polar angle, shows that the most essential contribution 
comes from |k| ~ |Xi2|~^, which, in view of inequality m2|Xi2| ^ 1, means that the 
suggested approximation is true. All this yields 



<^(;., |k|) ^ ( ^ ) ida.ida. 



/im| AvrVi, f, (Xi2-[X;xX'2]) 



-121 



247rT2 \ g ) J ' |X 

/ sinmz|Xi2| m\ \ 

The integral over the first term in the parentheses cannot be further simplified, but the sec- 
ond one can. Making use of the expansion Eq. ([4.7|) and extending the limits of integration 
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over z to ±00, in view of the fast convergence, one arrives at the expression for the damping 
rate: 



l-ln2 + M^(l-ll„2) + M^F. 



(4.17) 



with another geometric P-odd factor 
_X' • [X" X X'"] 



da- 



Xi2-[x;xx^] . 



A'iTLmz^i2 



sinmz X12 . 



(4.18) 



Note that for contours with a smooth dependence on a of the radius of curvature R , the 
first term in F2 is reduced to 



^daX.' ■ [n X n']/m|Li?2 



where n is the normal to the contour, and the integral over a is 27r times the so called 
twist [^] of the contour, while the second term, up to the factor 1/Lmz, looks similar to 
the writhing number (see them in Sec. [TVD , if it were not the oscillating factor in the 
integrand coming from the nontrivial distribution of Z-magnetic field. 



V. DAMPING OF THE BASIC CHARACTERISTICS OF Z-STRING 

Assume that the Z-string with the winding number no be formated. Then its energy 
is Eo ~ EvLuq. After that it goes to the state with the winding number n = 1 via either 
the fermion pair emission, after the time duration Tq ~ 2r~^lnno, or by some another 



mechanism mentioned in whose rate is unknown. In turn, the state with the single 
winding number decays only via the emission of the fermions. Since this state has the energy 
e^L and is the superposition of the stationary states whose excitation quantum numbers N 
are distributed according to the Poisson formula with the mean A'o = e^L/mz-, the number 
of emitted fermions is by a factor of In A^o ^ 1 grater than that in the case of the single Z° 
boson. Since Z-flux is directly related to the winding number, it decays with the rate r/2. 
The energy depends quadratically on Z-flux, so its damping rate is V. 

The energy and Z-fiux are nonzero for string configurations of arbitrary shape, including 
the straight strings. On the contrary, the helicity, 

d^k ik- [v(k,t) X v*(k,t)] _ 




(27r)3 (k^m^^ + 1)2 
\ 2 

^ j j daadab exp[-ik ■ (X^ - Xb)]nanb 



-g j J (2vr)3Vk2 + m|; 
x^k-[XlxX^]/k2 (5.1) 



2^], is nonzero only for the configurations of Z-strings which are not invariant under the 
space inversion. Indeed, the terms with a 7^ 6, after the momentum integration, give the 
linking number, 
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L[a, b] = — f d(7a f da, 



b- 



At^J ""J " X 



ab 



]26[] , of two contours, with the exponentially small corrections. The contribution of the 



typical term with a = b, after the momentum integration, reads 

hz{a = 6) oc W[a] - i- j da^ j da2^^^^^^^^ (l + mz|Xi2| + ^m||Xi2 

exp(-mz|Xi2|), (5.2) 
where X12 = Xa((Ti) — Xa(cr2) refers to the same contour a, and 

X12 ■ [x; X x'2] 



W\a\ = — (p da^ <h da- 



An J J X 
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2t\ is the writhing number of the contour a. The m^-dependent term in Eq. ( |5.2|) is 



evaluated with the help of the expansion ( [4.7| ) to give 

/rfaxi-[x>x:]. 



27rm| / 

In the case of sufficiently smooth contours the latter can be represented as ~T[a\/{mzRY-i 
where 

T[a] = — idaX! ■ [n x n'] 
27r J 

P^ l is the twist number of the contour a whose normal vector is n and the radius of 
curvature is R. So the twist contribution to the helicity is suppressed as {Rmz)~'^, and the 
resulting expression for the helicity can be written as [|,^,^ 



hz = (—] \j2nlW[a]+2j2n^n,L[a,b]\ . (5.3) 

\ 9 / \^ a a<b J 

Since hz is just Chern-Simons number of the Z-string field configuration, it characterizes 
possible processes with anomalous nonconservation of the baryon number 

The rate of the change in the helicity is evaluated semiclassically. To this end one 
should first take the expectation value of Eq. (|5.1|) in the quantum state discussed in 
Sec. |ITT[ Further differentiation with respect to time is then performed similar to the case 
of the damped classical oscillator, where the averaging over period is implied, giving the 
contribution due to the variable Z-flux to be = —Thz- The contribution coming from 

the slow classical contour motion is calculated with the help of the relation 

— ^ da[k X X'] exp(-ik ■ X) = i ^ rfak x (k x [X x X']) exp(-ik ■ X), 
which can be verified by a straightforward calculation. One finds 



g J J (2/1 ; yiv -r iiLz/ ab 

x exp(-zk ■ X„;,) 

^1 



27rm? 



9^ ab 



J^inaUb) j da a j> dabi±a - y^b)[^a X exp(-mz|X„b|). (5.4) 
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It is clear that the terms with a give exponentially small contribution. This is natural, 
since the analogous terms in the expression for hz give the contribution to the linking 
number L[a, 6] known to be the topological invariant. The contribution of the terms with 
a = 6 is calculated with the help of Eq. ( [4.7|) . In contrast to the case of helicity itself, their 
contribution to the time derivative of the latter is independent of niz- In total, one obtains 



The second term in the right hand side is of purely classical origin. Its inclusion is justified 
in the case when the decay by the fermionic pair emission is prohibited. 



The results obtained in this paper shed a new light on the issue of stability of the gauge 
vortex defects in that the latter can be considered as purely classical field configurations 
only under the definite circumstances. The winding number of the scalar field is intimately 
interrelated with the gauge field so that if the latter decays the former will do the same, 
provided the spectrum of excitations allowed to couple to the string background satisfies 
the specific threshold condition. The radial part / (modulus) of the scalar field is irrelevant 
in the London limit. The above threshold condition is satisfied for Z-strings coupled to 
fermions. The interaction with fermions of these extended objects is innate and will result 
in their decay even in the models that admits metastability at the classical level p!2| , p!3| . The 
suppression of this decay could be possible, in principle, at the nonzero fermionic density in 
surrounding matter. The main reason of instability of a sufficiently long Z-string in empty 
space is not the collapse (in the case of closed contour), but a much more fast process of 
the fermion emission. The numerical estimate for the critical length Lq is obtained from the 
condition 



Taking the velocity of the shrinkage v ~ 0.5, one gets Lqitlz ~ 2 x 10^. The loops with the 
length greater than 2000 Compton lengths of Z-boson will evaporate before their shrinkage. 




(5.5) 



VI. CONCLUSION 
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